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Abstract: This paper develops a European option pricing formula for frac- 
,-0 , tional market models. Although there exist option pricing results for a frac- 

O ■ tional Black-Scholes model, they are established without accounting for stochas- 

P_l ' tic volatility. In this paper, a fractional version of the Constant Elasticity of 

Variance (CEV) model is developed. European option pricing formula similar 
to that of the classical CEV model is obtained and a volatility skew pattern is 
revealed. 
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1. Introduction 

Black 0, 0] developed the so-called Constant Elasticity of Variance model for stock 
price processes that exhibit stochastic volatility. The CEV model is expressed in 
terms of a stochastic diffusion process with respect to a standard Brownian motion 



(1.1) dX t = fiX t dt + aXf /2 dB t , 

> . 

where < (3 < 2 is a constant. If /3 = 2, such a model degenerates to a geometric 
Brownian motion. This model is characterized by the dependence of the volatility 
rate, i.e., cV 3 / 2 on the stock price. When the stock price increases, the instanta- 
neous volatility rate decreases. This seems reasonable because the higher the stock 
[ price, the higher the equity market value, and thus the lower the proportion of 

liability, which results in a decrease in the risk of ruin. The volatility rate or the 
risk measure is thus decreased. Making use of methods proposed in an earlier liter- 
ature [|[ , Cox Q studied the CEV models and gave an option pricing formula that 
involves a noncentral x 2 distribution function. 

The classical CEV model (11.11) does not account for long-memory behavior, how- 
ever. There are some evidences showing that the financial market exhibits long- 
memory structures (see 0, HH). To encompass both long- memory and stochastic 
volatility, a possible model is to replace the Brownian motion in the stochastic 
■ diffusion equation by a fractional Brownian motion that exhibits a long-memory 

dependence structure (see [2, ljj 30j|)- 



Though fractional Brownian motion can be used to model long-memory, as 
pointed out by Rogers in 29], the fractional Brownian motion is not a semi- 
martingale and the stochastic integral with respect to it is not well-defined in the 
classical Ito's sense. A theory different from the Ito's one should be used to handle 



the fractional situation. One approach is white noise calculus (see [18|, l22|, |28(), 
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which was used in 14j, [l9( to construct stochastic integral with respect to the 



fractional Brownian motion. With the white noise approach, an extension to the 
Black-Scholes' stochastic differential equation is proposed to cope with long-memory 
phenomena (see [HI)- 

In this paper, the white noise calculus approach is adopted to construct a frac- 
tional CEV model and to derive the option pricing formula for European call option. 
Basic concepts in white noise calculus are briefly introduced in Section 2. The frac- 
tional Ito's lemma, which is fundamental to option pricing theory, is presented in 
Section 3. Section 4 explains under what circumstances is the Ito's lemma applica- 
ble. In Section 5, the fractional option pricing theory is introduced and the concept 
of self-financing strat egy , which is different from the traditional definition adopted 
by, e.g., Delbaen 13, HI, wm " further be discussed. Finally, the pricing formula for 



fractional CEV model is given in Sections 6. 



2. White noise calculus and stochastic integration 

In this section, the concept of stochastic integration with respect to fractional Brow- 
nian motion is introduced briefly. Important concepts are defined based on the white 
noise theory originated from fl7i | . who considered the sample path of a Brownian 
motion as a functional. Throughout this paper, notations used in [l, 14, 18, 2^ are 
adopted. 

Let S(R) be the Schwarz space. Take the dual fl = S'(R), equipped with the 
weak star topology, as the underlying sample space, i.e., to G f2 is a functional that 
maps a rapidly decreasing function /(•) g S(R) to a real number. Also, let B(Q) 
be the er-algebra generated by the weak star topology. Then according to Bochner- 
Minlos Theorem (see Appendix A of [HI), there exists a unique probability measure 
// on B(il), such that for any given / S S(R), the characteristic function of the 
random variable uj — > Lu(f) is given by 



n 

where 

||/|| 2 = / f(t)dt. 
Jr 

Let L 2 be the space of real- valued functions with finite square norm 1 1 • 1 1 , we have 
the triple S(R) C L 2 C S'(R). For any / S L 2 , we can always choose a sequence 
of /„ £ S(R) so that /„ — > / in L 2 , and uj(f) is defined as the \vm n ^, 00 Lo(f n ) in 
L 2 (p). ^ 

Consider the indicator function 

1, if < s < a, 
l(o,a)(«) = { -1, if a < s < 0, 
0, otherwise. 

It can be verified that for any two real numbers a and b, the random variables 
w (l(o, £!)(•)) and w(l(o ; 6)(-)) ar e jointly normal, mean zeros, and with covariance 
min(a, b). Define B(t) as w(l(o,t)("))> we can always find a continuous version of 
B(t), denoted by B, which is the standard Brownian motion. Roughly speaking, 
the probability space (ft, -B(f2), /x) can intuitively be considered as a space consisting 
of all sample paths of a Brownian motion. 
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Following the approach of we give the definition of fractional Brownian 
motion with Hurst parameter i < H < 1 in terms of white noise setting by using 
the fundamental operator Mjj, defined on the space L 2 , by 



M H f(x) 



ch J R f(t)\t-x\ H -idt,H < i < 1, 



where cr is a constant depending on the Hurst parameter H via 

c H = (2T(H \) cos(|(iJ - l)))- 1 [T(2H + 1) sin(TriJ)] * . 
Then, cj(Af//l( , a ) (')) an d uj(Mh l(o, &)(")) are jointly normal with covariance 

\{\ a r+\br-\a-br}. 

Again, we can find a continuous version of w(M^l( t )(-)) denoted by B (t), which 
is the fractional Brownian motion. 

We have the following Wiener-Ito Chaos Decomposition Theorem for a square 
integrable random variable on S'(R) (see Theorem 2.2.4 of |18|). 

Theorem 2.1. // F G L 2 (Q, B(Q), [i), then F(oj) has a unique representation 

F(o;) = ^c a JJ a (o;), 

ot 

where a is any finite integers sequence (a%, a>2, ■ ■ ■ ,a n ), c a are real coefficients and 
H a (uj) — h ai {uj{ei))h a2 (co^)) ■ ■ ■ h an (cj(e„)) h n {x) are Hermite polynomials and 
e n is an orthonormal set in S(R) which is defined as 

ei{t) = lyr^h^e-* 2 / 2 . 
Furthermore, the L 2 norm of the functional F(u>) is given by 

a 

with al = ai!«2! ■ ■ ■ a n l. 

Remark 2.1. (see [13, [S3l) The basis {H a {uj) : a} is orthogonal with respect to 
the inner product E(XY) in L 2 (Cl, £>(£!), /i). The variance of H a (uj) is al. Hq(uj) is 
taken as the constant 1. For a^O, the expectation of H a (cu) is 

E(H a (u)H (u))=0. 

As a result, the term Co is the expectation of the functional F(uj). 

Consider the functional Bf — aj(M^l[ 0j t] (.)), where t is a given constant. Using 
the dual property (see [l4j) of the Mh operator: for all rapidly decreasing functions 
/ and g, we have 

(f,M H g) = f f(t)M H g(t)dt= f g(t)M H f(t) dt = (M H /, g). 
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The function Muhorfi-) can be rewritten in the Fourier expansion form: 

oo 

M H l m (s)=J2(M H li ,t](-),e i (-))e i (s) 

oo 

= E( 1 [o,*](')' M ^ e ^-))e l (s) 

= ^2 <yj M H ei(u)dv)j e,(s). 
Since lo is linear, the functional can be written as 

oo 

w(M ff l [o , t] (-)) = X;(l[o,t](-).ei(0)w(ei) 

i=0 

= y2\ M H ei(u)du\ ui(ei). 

a 

In this example, when a = e(z) = {0..., 1, 0, ...}, i.e., one at position z, c a = 
J* MH£i{u)du, and c a ~ otherwise. It is tempting to write 

d 

-B? = 5^Mj r e i (t)w(e i ), 

i=0 

which is illegitimate in the traditional sense as the Brownian motion or the frac- 
tional Brownian motion is nowhere differentiable. With the chaos expansion form, 
differentiation and integration with respect to time t can be defined, but they may 
not always be square integrable. Such type of operation is called integration or 
differentiation in (S)* (see |18j). 

Definition 2.1. Let (S) be a subset of L 2 (il, B(il), /i) consisting of functionals 
with Wiener-Ito Chaos decomposition such that 

e I ^ n ( 2 ^ Qj 1 < °° 

a [ jeN J 

for all k < 00 and that (5)* consists of all expansions, not necessarily belonging to 
L 2 (tt,B{tt),Li), such that 

a [ jeN J 

for some q < 00, then, the spaces (S) and (S)* are called the Hida test function 
space and the Hida distribution space respectively. 

The derivative of the fractional Brownian motion, or the white noise is defined 

by 



W H (t) = Y j M H e i {t)u(e i ). 
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It can be shown that the sum W H (t) G (S)* (see [14|,|3lJ). The importance of Hida 
test function space and distribution space is their closedness of Wick multiplication. 
Wick product is an operator acting on two functionals F(co) and G(lo). 

Definition 2.2. The Wick's product for two functionals having Wiener-Ito Chaos 
Decomposition 



and 



is defined as 







F(ui) o G(w) = c a bpH a+l3 (uj). 



Addition of indexes refers to pairwise addition. 

The closeness of Wick's product is shown in the following theorem (see Corollary 
2.2 and Remark 2.8 of [Hj]). 

Theorem 2.2. Wick multiplication is closed in (S) and (S)* . 

It is reasonable to define the stochastic integration of a functional Z t (cu) with 
respect to the fractional Brownian motion as the integration with respect to time 
t of the Wick's product between Z t (uj) and W^{u>). Under the Wiener Chaos de- 
composition framework, if the decomposition exists, the functional Z t (io) o (to) 
can be written as 

Qi 

It is natural to think that the integration is 



E 



(i 



c a (s)ds \ H a (u), 



by assuming that summation and integration are interchangeable. If the integration 
with respect to time is a path-wise classical Riemann integral, it is not clear that 
summation and integration are interchangeable. Such difficulties can be finessed 
by introducing new definitions for integration with respect to time and integration 
with respect to the fractional Brownian motion as follows (see Definitions 2.3 and 
2.4 respectively). 

Definition 2.3. (a) Elements in (S)* as an operator: Let F(u) — c a H a (u) £ 
(S)* and f(ui) — J2 a b a H a (^) & (S), then F can be regarded as an operation on / 

(FJ) =J2b a c a a\. 

ot 

(b) Time integration: If F t {ui) = J2 a c a{t)F[a{co) are elements in (S)* for all 
positive real number t, and that (Ft, f) are integrable with respect to t for all 
/ G (S), then the integral f R F t (ui)dt is defined as the unique element in (S)* , I(w) 
such that 

(I(uj)J) = f (F t (u),f)dt 

for all / G (S). 
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Remark 2.2. It can be shown that the quantity (F, f) in part (a) exists and is finite 
under the condition given in the definition. It can be regarded as the expectation of 
the product between F(io) and f(u>) when F(u>) £ L 2 (fl, B(fl), fi) and f(w) £ (S). 

Definition 2.4. Let Z(t) = J^a^i^Hai^) € (£)* for any given t, then, the 
Wick's integral of Z(t) is defined as 

f Z(t)odB H (t)= [ {Z(t)oW H (t)}dt, 

JR JR 

when Z{t)o W (t) is integrable with respect to time t in the sense as in Definition 
2.3. 

The following theorem asserts that integration and summation are interchange- 
able, see Lemma 2.5.6 of (lif . 

Theorem 2.3. Let Z(t) : R — + (5)* , with Wiener Chaos decomposition 

y^^c a (t)H a (u) 

a 

such that 

J2 a] -{[ c a (t)dt} 2 Y[(2j)-^ <co 

for some q < oo, then Z{t) is time-integrable, also, integration and summation are 
interchangeable, i.e. 

I Z{t)dt = Y\{( c Q (t)dt}H a {uj). 
JR a JR 

3. Fractional Ito's lemma 

Several approaches of extending classical Ito's lemma to incor por ate the fractional 
Brownian motion were discussed in the literature, e.g., foL Il2j|. The settings in 
these papers are different and various conditions are required to ensure that the 
stochastic integrals appear in the fractional Ito lemma exist. Bender in \1\ provided a 
simpler version of fractional Ito's lemma based on the white noise setting introduced 
in the preceding section. Here, we restate Bender's theorem in Theorem 3.1 and 
give a generalized result in Theorem 3.2. 

Theorem 3.1. Consider the stochastic process 

Y t = f h{t)dB» =u[M H (h(-)l m (■))], 
Jo 

where h(t) is a continuous function in [0, T] and H > h. Let g(t,y) be a two- 
dimensional function differentiable with respect to t and is twice differentiable with 
respect to y. Also, there exists constants C'i > and 
Ai < (2T H sup se r T i h(s))~ 2 so that 



max{|fif|, Ifftl, \g y \, \g yy \} < Cie Al 
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Let k(t) — H(2H — l)/i(t) J Q h(s)\s — t\ 2H 2 ds. Then, we have the following frac- 
tional Ito's lemma: 

g(T,Y T ) = 5 (0,0) + J ^g(t,Y t )dt + J h(t)^-g(t,Y t ) o dBf 

+ J k{t)—g(t,Y t )dt. 

Remark 3.1. The integrals above are well defined as the condition that the inte- 
grands and integrals both belong to (L 2 ) is ensured by the given assumption. Since 
(I/ 2 ) C (S)*, the integrals are defined as in Section 2. 

This theorem gives the differential form of g (t, Y t ) when the underlying stochastic 
process Y t is an integrals of a deterministic function h t . The following generalization 
takes the underlying stochastic process to be X t = g(t, Y t ) and gives the differential 
form for P(t,X t ), where P(t,x) is a two-dimensional real- valued function. Before 
introducing our results, let us illustrate some ideas through an example (see (Hf). 

Consider the two-dimensional function, 

g(t, y) = exp^i - \o 2 t 2H + ay), 
where fj, and a are two positive constants and the underlying stochastic process is 

Y t = [ dB*, 
Jo 

i.e., h(t) — 1. Clearly, for any given value of T, the functions g 7 g t , g y and g yy are 
all continuous in the closed interval [0, T] and hence, the conditions in the theorem 
are fulfilled. Applying the lemma, we have 

k(t)=H(2H-l) f \t- s\ 2H - 2 ds = Ht 2H -\, 
Jo 

and 

dX t = (fi - ^Ht^-^git, Y t ) dt + crg{t, Y t ) o dBf 
+k(t)a 2 g(t,X t )dt 
= t ig(t, Y t )dt + ag(t, Y t ) o dBf 
= [iX t dt + o-X t odBf. 

The next question is whether there exists an Ito's lemma that further expresses 
P(t,X t ) in terms of integrals involving fiX and o~X, but not Y and g explicitly. It 
is reasonable to expect that 

dP(t, X t ) = P t (t, X t ) dt + tiY t P x (t, X t ) dt + o-Y t P x (t, Y t ) o dBf + 
<j 2 Ht 2H - 1 X 2 P xx (t,X t )dt. 

This result can be verified by the following theorem. 

Theorem 3.2. Using the same notations and assumptions in Theorem 3.1, further 
assume that the differential of the stochastic process X t = g(t,Y t ) can be, according 
to Theorem 3.1, written as 

dX t = n{t, X t )dt + a(t, X t ) o dBf. 
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Also, there exists constants C2 > and A2 < (2T H sup se r T i /i(s)) 2 so f/iai the 
composite function P o g = P(i, g(i, y)) satisfies 

max{|(Po 5 )|,|(Fo 5 ) t |,|(Po 5 ),|,|(Fo 5 ) w |}<C 2 e A ^ 2 . 
Let C(i) = k(t)h~ 2 (t). Then the fractional Ito's lemma is given by 

P(T, X T ) = P(0, X ) + £ ^P(t, X t ) dt + £ n(t, Yt)-^P(t, X t ) dt 

+ J a ^ Y t)^ P ^ X t)odB^ + J C(t)a 2 (t,Y t )^P(t,X t )dt. 

Proof. This theorem can be verified by applying Theorem 3.1 to f(t,g(t,Y t )). □ 

To relate this result with previous works on fractional stochastic calculus and 
option pricing theory, such as 12 , 2(| , consider the last term in the right-hand side 
of the Ito's formula. This second order correction term can be considered as the 
product of three quantities: P xx (t,X), a(t,X) and C(t)a(t, X). Comparing with 
the result of [13], the quantity C(t)a(t, X) corresponds to Dfx t . This quantity 
is known as the Malliavin derivative of X t . For details of Malliavin calculus, see 
2oL 24, 23]. Under our assumptions that X t has the form of g(t, Y t ), this quantity 
depends only on the current time t and the value of X at time t. But this needs 
not be the situation for a general X t governed by a fractional stochastic differential 
equation. In general this quantity may depend on the entire path of X t , not only 
on the value at one point. This may introduce further complications when working 
with the differential. 



4. The fractional CEV model 

Here we construct a fractional version of the Constant Elasticity of Variance model 

by means of the fractional Ito's lemma (Theorem 3.1). As discussed in the intro- 

3/1 

duction, constant elasticity is characterized by the the volatility term o~X t in 
a stochastic diffusion equation. In order to handle long-memory, we replace the 
Wiener process by a fractional Brownian motion. The fractional diffusion equation 
is then defined as 

dX t = fx(t, X t ) dt + crXf 12 o dBf , 

where 0</3<2. If_ff = i and /j,(t,X t ) = [iX t , this is the classical CEV model. 
In this situation, the Wick integral is equivalent to the Ito's integral (see fljj ) 
and hence, the classical Ito's lemma can be applied to any stochastic process of 
the form Yt = P(t,Xt). When long-memory is considered, the Ito's lemma will 
involve the Malliavin derivative, which is in general path dependent and difficult 
to handle. When the integration is defined in the white noise sense (Section 2), 
the integrand is assumed to belong to (£)* at every time t, and the integral is a 
random variable in (S)* . The elements in (S)*, which are merely formal expansions, 
may not correspond to real values for each lo and the term P(t,X t ) may not be 
well-defined in general. In order to overcome such difficulties, we need to choose a 
suitable /j,(t, X t ). 

Assume that X t can be written explicitly in terms of time t and a stochastic 
integral process Y t = J*h(s)dB^, i.e. X t = g(t,Y t ). From Theorem 3.1, the dif- 
ferential of X t can be decomposed into two parts, the drift term and the volatility 
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term. In order to keep the elasticity constant, the function g(t, y) must be chosen 
as 

9(t,y) = L(l-^)[h(t)]- 1 y + f(t) 

where f(t) is an arbitrarily chosen function of t. This g(t, y) is in fact the solution 
to the differential equation 

h(t)^-g(t,y) = <j{g(t,y)}^. 
dy 

The next question is how to choose fit). The answer is given by the following 
theorem. 

Theorem 4.1. Assume that h(t) is a strictly positive function and g(t,y) is the 
solution to the differential equation 



(4.1) 



h(t)-^-g(t,y) = a(t,g(t,y)). 



The general solution to this equation involves an arbitrary function fit). 
Let r)(t) and (p(t) be two functions determined by the integral equations, 



h(t) = e 



and 



ao + / h(s)<p{s) ds 



fit) = [hit)]- 1 
where ao is a constant so that g(0,0) = X$, then, 

X t =g(t,Y t ) 

yields the volatility o~it,X t ) in the fractional Ito's lemma and the drift term is given 
by a two dimensional function 



ax JO o\p,X) 



ait, x) 



dx 



ait, x) 



Hit, x) = ait, x) 
Proof. Let 

then git, y) given by 

ait,git,y)) = [hit)}- 1 y + fit) 
satisfies Equation (|4.1[) . After some manipulations, we have 



1 da ' 
^— dx 



a 2 it,x) dt 



dt 



. . . da , . . . . „ da 

v{t)+C{t) — {t,g) + r,{t)ait,g)- — 



which does not involve y explicitly. The proof is completed by comparing this with 
fractional Ito's lemma (Theorem 3.1). □ 



158 



N. H. Chan and C. T. Ng 



Remark 4.1. Note that the function C(t) depends on the choice of r){t) and by 
definition (in Theorem 3.2) must be strictly positive. 

In the CEV case, when the drift term has the form of 

(4.2) x) = J A-x + Mt)x^ + °^x^C{t), 

the solution to the stochastic diffusion equation 

dX t = fj,(t, X t ) dt + crXf /2 o dBf 

is well defined and is given by Theorem 4.1. Now choose suitable time dependent 
functions r](t) and (p(t). Consider the three quantities of /i(t,X t ) in (|4.2p . The last 
one is strictly positive and cannot be eliminated. When we choose 77(f) = (1 — ^)fi 
and ip(t) = 0, the second term vanishes and the first term becomes /J.X t . The drift 
term becomes 

»{t,x t ) = nX t + ( ^-c{t)xt 1 . 

The CEV stochastic differential equation is thus 

X t = X Q + nX s ds + ^C{s)XP- 1 ds + J aX^ 2 o dBf X . 

5. Fractional option pricing theory 

By using the generalized Ito's lemma (Theorem 3.2), the differential of P(t, X t ) can 
be decomposed into two terms, the drift one and the volatility one and both terms 
involve only current time t and X t , i.e., they are path- independent. The foundation 
of the Black-Scholes' option pricing theory is constructing a self-financing strategy, 
which makes use of stocks and bonds to hedge an option. The definition for self- 
financing strategy in continuous-time model depends on how the stochastic integrals 
are defined. As the fractional stochastic integrals arc defined in a different manner, 
a new definition for self-financing strategy is required. One approach is adopted 
by 14, [n|, which defines self-financing strategy under the geometric fractional 



Brownian motion. Here, this approach is extended to a more general situation. 
Let Xt be the stock price process and lit be the bond value and they are governed 

by 



X t = [ n(s,X)ds+ [ a{s,X)odB^, 
Jo Jo 



(5.1) II t = / rllds 







Definition 5.1. (see Section 5 of |14[) A trading strategy consists of a quantity 
(u t ,Vt) of bonds and stocks is called self-financing if the infinitesimal change in the 
portfolio value at time t is given by 

dZ t = d(u t Ut + v t X t ) 

= rn t ut dt + fi(t, X t )v t dt + [a(t, X t )v t ] o dBf + rfA, 



where dA. is an infinitesimal dividend payment term. 
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Below, using the above definition, a fractional Black-Scholes equation is derived. 

Theorem 5.1. Suppose that the market consists of two securities, a risk-free bond 
and a stock. Here, the stock provides dividend continuously with rate S. Assume that 
the stock price process X t = g{t, Y t ) is defined as in Section 3 which also satisfies 
the equations i5. Then the price of a derivative on the stock price with a bounded 
payoff f(X(T)) is given by P(t,X), where P(t,X) solves the PDE: 

BP B 2 P BP 

(5-2) ^ + o\t, X)C(t)— + (r - S)X — rP = 0, 

with boundary condition P(T,X) = f(X) given that P o g satisfies the conditions 
in Theorem 3.2. 

Proof. A proof parallel to the fractional Black-Scholes equation in (2(| is given here. 
Consider a solution P(t,X) to equation (|5.2| . Applying fractional Ito's Lemma 
Theorem 3.2, 

BP dP d 2 P BP 

dP(t, X t ) = gx»(t, X) + Q^cj\t, X)C(t)} dt + —a(t, X) o dBf. 

Form a trading strategy by dynamically adjusting a portfolio consisting a varying 
quantity v[t) of stocks and u(t) of bonds. By choosing 

, , dP 

»(" = sx 

(5.3) „( () = j L ( P-x|f), 

then the portfolio value at time t is Pt and 

ru(t)n t dt + v(t)n(t, X t ) dt + [v(t)a(t, X t )\ dB? + Sv(t)X t dt 

BP BP BP BP 

= [rP rX—} dt + n— dt + (a—) o dB? + — SX dt 

BP BP B 2 P 

=^- d x 5xdt+ 9x-y^ c ^ dt 

BP BP BP 

+ dt + (<r%-) o dB? + ^—SXdt 
f 3X v OX' 1 dX 

BP BP B 2 P BP 

= [ -bJ + bx^ x) + Bx^ a2[t ' x)cm dt + dx a{t > x) ° dB " 

= dP(t, X t ) 

= d{u t n t + VtXt). 

By Definition 5.1, (u(t),v(t)) is a self-financing strategy. It can be shown that such 
strategy hedges the derivative. The portfolio value at time t is given by u(t)H t + 
v(t)X t and it is equal to P(t,X t ). At time of maturity, the portfolio value is just 
P(T,Xt). By assumption, the function P(t,X) satisfies the boundary condition, 
so P(T,X T ) = f{X T ). Therefore (u(t),v(t)) hedges the derivative and P(t,X) is 
the option price. □ 

6. Pricing an European call option under CEV models 

Putting a 2 (t, X) = a 2 XP, the Black-Scholes PDE (Theorem 5.1) of the CEV model 
is now given by 
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Putting Y = X 2 - f3 (see d) and P(t, X) = e rt Q(t, Y), this equation becomes 

^ + [bY + cC(t))^+aC(t)Y^=0, 

where a = a 2 (2 - fj) 2 , b = (r - S)(2 - fj) and c = a 2 (2 - 0){l - fj). The boundary 
condition is 

Q(T,Y) = e" rT max(y^,0). 

The approach of Cox and Ross Q that made use of Feller's result ( TH, can 
be adopted. First the solution for this equation at (to,Y to ) is the expectation of 
Q(T, Y T ) under the SDE 



(6.1) dY = [bY + cC(t)] dt + y/2aC(t)Y dB u 

with Y(to) = Y to (see 0). To solve this SDE, we follow Feller's arguments. First, 
a useful result of Kolmogorov is stated below (see Equation (167) of [ll[). 

Theorem 6.1. The probability density function of a diffusion process X t driven by 
standard Brownian motion 

dX t = fi(t, X t ) dt + a(t, X t )dB t 

is given by the PDE 

A 



u t = [-a 2 (t,X)u(t,X)] xx - [»(t,X)u(t,X)] 



x- 



L 2 

In our case, because of (I6.1[) . the Kolmogorov's equation becomes 

u t = [aC(t)Yu(t, Y)] YY - [(bY + cC(t))u(t, Y)] Y . 

The European call option pricing formula can be obtained by solving the above 
PDE. 

Theorem 6.2. Under the fractional CEV model introduced in Section 4, the price 
of an European call option with strike price K , mature at T at current time to is 
given by 

P(t , Xo) = e-^- f »»l„ £ ± e -^(^YG(r + 1 + * ^-) 

r! ajT 2 — fj ae ul 7t 

- Ke-^~^ J2 w-T i S e-^(—) r+ ^G(r + 1, -^), 
f^ Q r(r + 1 - 2^a) «7t ae bI lT 

where x = e" w °Y(t Q ), 

Jt = f e- bT C(r)dT 

and 

1 f 00 

r(a) ,/„ 

Proof. Assume that the Laplace Transform of u(t, Y) with respect to Y exists and 
equals to uj(t, s). Since the value of Y at time to is given, Y to is deterministic and 
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thus u(to,Y) = S(Y — Yt ), the Dirac function. Also, ui(to,s) = e sY " = tt(s) and 
the equation becomes the boundary value problem 

(6.2) u t + s(aC{t)s - b)uj s = -cC(t)su + ip(t) 

(6.3) uj(to, S ) = e- sY °. 

In equation (|6.2[) , ^(-) is called the flux of u at the origin (see [IB|), which is to 
be determined later. Now, we find the characteristic curve of the first order PDE 
(|6.2p . The characteristic curve is given by 

^ - aC(t)s 2 - bs. 

The solution to this equation is 

s = e- H [Ci - a lt ]-\ 
On this curve, u>(t, s(t)) satisfies 

-u{t, s(t)) = i/,(t) - cC{t)s(t)uj(t, s(t)). 

Solving, we have 

*(*)) = [C 1 -a lt ] c t a [C 2 + f V(r)|Ci - a lT \- c ' a dr]. 

For any given point [t\, s\), the characteristic curve with 

Ci = 07tj H ~T7~~ = C(ti,si) 

Sie otl 

passes through (ti,si). Also, this curve passes through the point (to, C^ 1 e^ bt °). 
This yields the value of C 2 , 

C 2 = [C 1 (t 1 ,s 1 )]- c / a u(to,e- bt °Ci 1 (t 1 ,s 1 )). 

The Laplace transform of u(t,Y) at point (ti, Si) is thus given by 

u){t U 8x) = (s 1 e bt ^ a [(C 1 (t 1 ,s 1 ))- c / a cj(t ,e- bt °C^(t 1 ,s 1 )) 

+ ( 1 ip{r)\a{ ltl - lT ) + J—\-e/»dT] 



6(ti-to) 

[W* 1 ^ + 1]- C /V(^i- — ) 

\siae b *i7i 1 + V 

+ [ 1 { Sl ae bu ( 7tl - 7t ) + l]- c/ >(r) dr. 



(6.4) 



J to 

Following the argument of [ljj] . when «(t, 0) is finite and c<0or0<c<a, 

lim(sae b *7 t + lWi, s) = 0, 

S— > 

then is given by the integral equation 

7r( — ) + [\{t){— l^) c / a dT = 0. 

07t Jto 7t - 7r 
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To solve this equation, applying the substitutions z = — and C = — , 



(6.5) 



.9(C)(C - *r c/o rfC = A — ) = A ) 



where g(() = -^( T )C /a %- 
The solution of (16.51) is 



: exp( ), 

a 



3(0 = 1^1)$^ «p(-^0. 



^(r)=9(C)r /o f 



1 , 1 w a; . -o+q , x dj T 
(— )( ) ° ex P( )" 



r(l - |) 7 r a7 T a7 T dr 

Substituting this result into (|6.4[) . after simplification, we get 



w(t, a) = cxp( r^-)( ^ -) C / Q r(l - -; - % -). 

The next step is to perform an inverse Laplace transform with respect to s. To this 
end, let 

A= — , 

2 = sae bt j t + 1- 

One can verify that equation (6.5) in [16| | is still valid after these substitutions. The 
quantity w(t, s) can now be rewritten as 



-A A 1 



Jo 



r (!-f) 

Using the fact that Laplace Transform of Iq(z(AtY) 1 / 2 ) is e^z" 1 , we have 

u(t,Y) = ( * )( ^ ) ^exp{-^±5^}J 1 _ f [A (e -^y)V 2]) 
ae 0I 7t r ae ot 7 f a a7 t 

where /a(-) is the first type Bessel function with order A, which is defined as 

- ( . /2) 2*+A 

^ /c!L(fc + 1 + A) ' 

This density function is then used to find the solution of P at (to, Xq) by means of 
the identity, 

P(t ,X ) = e rt "Q(t Q ,X ) = e-^-^^M^ - K,0)}. 
After direct calculations, we have 

/•OO 

P(t , X ) = e -< T -^ / (y*h - K)u(T, y) dy 



IK 2 

1 __«_, a: 1 X 2 "' 3 



r(T _ t0) £ e (,- 5 )T^ J_ ^ ( _E_)r G(r + 1 + - J-, -ii^) 

^ — ' r! a7T ^ — p ae 7t 



; e 

r=0 



^T(r+ 1-5^3) K a lT ae bT lT 
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Putting x = e 



bt y2-0 



o > 



x 



2^ = e'^Xo = e-^-^X, 



So the option pricing formula is 



P(t Q ,X Q )=e- 5 ^X Q Y, 




OO 



i> + i-^) e 



e -^7(^)''+2 J -?G( r+ l ) 



)• 



r=0 



ae" 1 7t 



□ 



This formula is similar to the classical one, which is obtained by replacing the 
term 77- by the term ^ b -(e~ bt ° — e~ bT ). As these two terms do not depend on the 
strike price, the implied volatility pattern is the same as the classical CEV model. 
Consequently, the fractional CEV model also accounts for the volatility skewness 
observed in practice. 
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